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Abstract 

Henry's law constants of the solutes methane, nitrogen, oxygen and carbon dioxide in the 
solvent ethanol are predicted by molecular simulation. The molecular models for the so- 
lutes are taken from previous work. For the solvent ethanol, a new rigid anisotropic united 
atom molecular model based on Lennard-Jones and Coulombic interactions is developed. 
It is adjusted to experimental pure component saturated liquid density and vapor pressure 
data. Henry's law constants are calculated by evaluating the infinite dilution residual chem- 
ical potentials of the solutes from 273 to 498 K with Widom's test particle insertion. The 
prediction of Henry's Law constants without the use of binary experimental data on the ba- 
sis of the Lorentz-Berthelot combining rule agree well with experimental data, deviations 
are 20 %, except for carbon dioxide for which deviations of 70 % are reached. Quantitative 
agreement is achieved by using the modified Lorentz-Berthelot combining rule which is 
adjusted to one experimental mixture data point. 
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1 Introduction 



The goal of the 2nd Industrial Fluid Properties Simulation Challenge 2004 organized by the 
American Institute of Chemical Engineers (AIChE) is "to assess the current abilities and in- 
abilities in the prediction of physical properties applying force fields and molecular techniques" 
One problem proposed in this Challenge is the prediction of Henry's law constants for the 
solutes methane, nitrogen, oxygen and carbon dioxide in the solvent ethanol at the tempera- 
tures 323 and 373 K. This problem is tackled here using multi-center Lennard- Jones, polar and 
electrostatic interaction potentials together with molecular simulation. 

For modeling the solutes, the symmetric two-center Lennard- Jones plus pointpolarity po- 
tential is used. The parameters of this model were adjusted for 80 pure fluids to experimental 
saturated liquid density and vapor pressure data in previous work of our group (21131. 

Many molecular models for ethanol are available in the literature, some of them consider all 
atoms explicitly, some of them use the united atom approach for methyl and methylene groups. 
Most of the existing models account for internal degrees of freedom, but only few explicitly for 
polarizability effects. An overview of ethanol models from the literature is given in section [2?n 

Among the available molecular models for ethanol, the probably most appropriate for pre- 
dicting vapor-liquid coexistence properties is the transferable potential for phase equilibria- 
united atom (TraPPE-UA) of Chen et al. JH. That model accounts for internal degrees of free- 
dom. In the present work, a new simple molecular model for ethanol was developed with the 
aim to yield at least as accurate results but with distinctly less computational and programming 
effort. The new rigid effective ethanol model is of the united-atom type and based on Lennard- 
Jones and Coulombic interactions. The united atom approach and the neglect of internal degrees 
of freedom saves computation time which is desirable for applications of force field methods in 
industrially relevant problems. 
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An infinitely diluted solute in a solvent only experiences solute-solvent interactions. The re- 
sults for Henry's law constants from molecular simulations therefore depend on the model used 
for describing the unlike solute-solvent interaction. Consequently, predictions of the Henry's 
law constant from pure component data alone are a hard test for every molecular model. The 
unlike interactions are usually determined from the like interactions which are known from pure 
component models by combining rules like that of Lorentz-Berthelot. These combining rules 
usually only have a weak theoretical basis so that it is an open question how useful they are for 
predictions of Henry's law constants from pure component data alone. On that background it 
makes sense to also study modified combining rules which allow an adjustment of the model 
for the unlike interaction to experimental mixture data. Experience shows that if such an adjust- 
ment is carried out even for only one single mixture data point, excellent predictions of mixture 
properties can be achieved over a wide range of states Il51-[|8l. 

Henry's law constants were predicted by applying the molecular models which were param- 
eterized exclusively to pure substance vapor-liquid equilibria and by using the Lorentz-Berthelot 
combining rule for the unlike Lennard- Jones interactions. Also some other combining rules are 
briefly discussed. Additionally, the modified Lorentz-Berthelot combining rule with one state 
independent binary adjustable parameter is investigated. 

2 Molecular Models 

2.1 Ethanol models from the literature 

In this section, an overview of molecular ethanol models from the literature based on the 
Lennard- Jones potential and point charges is given. The models are assessed regarding their 
description of vapor-liquid equilibria where such information is available in the literature. In 
1981, Jorgensen |l9l presented a rigid ethanol model based on transferable intermolecular poten- 
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tial functions (TIPS). TIPS consider the methyl and methylene groups as single Lennard- Jones 
sites, centered on the carbon nuclei and a single Lennard- Jones site centered on the oxygen 
nucleus accounting for the dispersion and repulsion of the hydroxyl group. Point charges are 
located on the hydroxyl hydrogen and the bonded oxygen nucleus. The charge located on the 
methylene group is chosen to achieve overall neutrality of the molecule. That model was devel- 
oped to yield reasonable structural and energetic results. Subsequently, Jorgensen investigated 
the effect of internal rotation in the ethanol molecule about the carbon-oxygen bond IfTOll . re- 
garding the structure and hydrogen bonding and compared the results to X-ray data for the solid 
and to ab initio molecular orbital calculations. 



Improvements of the Lennard- Jones parameters of the united alcyl [flTll and hydroxyl groups, 
the point charges, the internal rotational potential and the geometry, were achieved by fitting di- 
rectly to experimental thermodynamic and structural data, as well as quantum chemical molec- 
ular mechanics calculations. This yielded the united-atom optimized potential model for liquid 
simulations (OPLS-UA) lfT2l . Vapor-liquid equilibria with the OPLS-UA were calculated by 
van Leeuwen (T3\ . A comparision of the these simulation results with experimental data shows 
a mean unsigned error [IJ of about 30 % in vapor pressure, 6 % in saturated liquid density and 
7 % in heat of vaporization. Stubbs et al. [[T4| suggested a modification of the Lennard- Jones 
parameters of the OPLS-UA methyl group. Petravic and Delhommelle IfTSl termed this modi- 
fied model S2 and compared it to the OPLS-UA regarding liquid density, hydrogen bonding and 
diffusion. From that study, it was concluded that the S2 model improves the prediction of liquid 
density. However, the S2 model overestimates the diffusion coefficient by approximately 15 %, 
whereas the OPLS-UA gives better results. Also the influence of the internal degrees of free- 
dom on transport quantities has been studied by Petravic and Delhommelle. They showed that 



1 1 ^ Iz^-z^"'"! ■ exn 

mean unsigned error; ^ exp — , where M is the number of simulated data, z-™ is the simulated and Z; 

!= 1 

the experimental quantity of interest. 
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the density is affected only slightly by the internal degrees of freedom, in contrast to diffusion 
coefficients. 

More detailed potential models for ethanol consider all atoms explicitly. Miiller-Plathe (TE] 
designed such an all atom model including internal rotation about the carbon-oxygen bond to 
give reasonable bulk properties in mixtures with water. Jorgensen et al. [17J published an OPLS 
all atom potential model (OPLS-AA) considering additionally stretching, bending and internal 
rotation about the carbon-carbon and carbon-oxygen bond. A similar molecular model which 
uses the same approach as OPLS-AA but with a different parameter set is the AMBER94 force 
field [fT8l[T9ll . Panagiotopoulos EOl concluded that the OPLS models are not appropiate predict- 
ing vapor-liquid coexistence since they were developed for energetic and structural properties 
at ambient conditions. 

Optimized models for vapor- liquid equilibria of linear and branched alkanes were developed 
by Martin and Siepmann [|2T1 122|| . i.e. transferable potentials for phase equilibria-united atom 
(TraPPE-UA). The TraPPE-UA was extended to alcohols ['4| which gives to our knowledge the 
most accurate existing ethanol model for vapor-liquid equilibria. It was found that this model 
yields vapor pressures, saturated liquid densities and heats of vaporization with mean unsigned 
errors of about 8, 1 and 1 %, respectively. 

Since polarazible molecular models are computationally very expensive, only few ethanol 
models of that type are available. The polaraziable intermolecular potential function (PIPE) 
proposed by Gao et al. [|23ll is based on the reparamerization of the OPLS-AA plus isotropic 
point polarization on each atom. That model was developed to yield better results for heat of 
vaporization and density. Another polarizable model is that from Gonzalez et al. Il24l . It is 
based on the OPLS-UA and consists of only one isotropic point polarizability on the (pseudo) 
oxygen atom. Delhommelle et al. [25] showed for pentane and hydrogen sulfide that polarizable 
models can give better predictions of mixture vapor-liquid equilibria. 
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2.2 New ethanol model 



A new ethanol model was developed with the aim to give accurate vapor-liquid equilibria with 
low computational and programming effort. It neglects the internal degrees of freedom and uses 
nuclei off-center Lennard- Jones united atoms for the methyl, methylene and hydroxyl group, 
accounting for repulsion and dispersion. Point charges are located on the methylene and hy- 
droxyl Lennard- Jones centers, as well as on the nucleus position of the hydroxyl hydrogen. The 
Coulombic interactions account for both polarity and hydrogen bonding. 

The molecular interaction and geometry of the resulting model is described in the following. 
The potential energy Uij between two ethanol molecules i and j is given by 

4 4 
a=lb=l 

where a is the site index of molecule i and b the site index of molecule j, respectively. The site- 
site distance between molecules / and j is denoted rijab- Oab, ^ab are the Lennard- Jones size and 
energy parameters, qta and qjt, are the point charges located at the sites a and b on the molecules 
i and j, respectively. Finally, Eq is the permittivity of vacuum. The interaction between unlike 
Lennard- Jones sites of two ethanol molecules is defined by the Lorentz-Berthelot combining 
rules 



<3ab 

^ ijab 
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(1) 



<3aa+<3bb 
<^ab = 2 ' 

^ab = V^aa^bb- (3) 



The geometry of the present ethanol model and the potential parameters are given in Figure \T\ 
and Table m 
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In the following, a few remarks on the model development are given. Jorgensen's inves- 
tigations on the potential of the internal rotation about the carbon-oxygen bond in an ethanol 
monomer [flOil yield three energetic extremal values in the corresponding configurations: the 
absolute minimum is reached in the trans configuration (dihedral angle (|) = 180°), the 0.64 
kcal/mol higher local minimum in the gauche configuration ((|) ^ 60°) and the 2.06 kcal/mol 
(above trans) higher absolute maximium in the cis configuration ((|) = 0°). For the present rigid 
ethanol model, the energetically most likely trans configuration has been chosen. The nuclei 
positions of all ethanol atoms were computed using the quantum chemistry software pack- 
age GAMESS (US) The basis set 6-3 IG and the Hartree-Fock method were applied for 
geometry optimization. The distances between the carbon-carbon, carbon-oxygen and oxygen- 
hydroxyl hydrogen nuclei positions obtained from quantum chemistry calculations are 1 .5 1342, 
1.47986 and 0.95053 A, respectively. The CCO and COH angles are 110.760 and 117.280°. 
Figure [T] includes also these nuclei positions. 

Starting from these positions, for the methyl and methylene group the AUA4 parameters of 
Ungerer et al. [|27ll were applied. The AUA4 parameters were optimized by Ungerer et al. for 
vapor- liquid equilibria of linear alkanes. Following the approach of Ungerer et al. [|27ll , a small 
offset of the Lennard- Jones hydroxyl center (0.1010 A) in direction to the hydroxyl hydrogen 
nucleus was allowed for optimization. The five parameters of the hydroxyl group were fitted to 
yield optimal saturated liquid densities and vapor pressures. The parameters are the two point 
charges, the Lennard- Jones size and energy parameters as well as the offset of the hydroxyl 
Lennard- Jones center. These parameters were chosen since they model the strongly interacting 
hydrogen bond and since the methyl and methylene parameters were already optimized by 
Ungerer et al. fTl]. The point charge on the methylene center was set in such a manner to yield 
overall neutrality of the molecular model. 

The proposed set of charges in the present work yields a resulting dipole moment of 2.47 D. 
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Quantum chemistry calculations (M0ller-Plesset level 2 with basis set aug-cc-pVTZ) in a cavity 
(COSMO with a dielectric constant £]■ = 25) and in the vacuum yield dipole moments of 1.99 
and 2.10 D, respectively. Note that the higher dipole moment of the molecular model than those 
from quantum chemistry calculations is not unusual for modeling hydrogen bonding molecules 
with point charges. 

Pure substance vapor-liquid equilibria for the optimization of the ethanol model were cal- 
culated using the NpT+test particle method proposed by MoUer and Fischer [|28l - [l30l . To 
determine the residual chemical potential in the liquid phase with high accuracy, the Monte 
Carlo based gradual insertion method OTI - OSlI was used. For the residual chemical potential 
in the vapor phase, standard Monte Carlo simulations with Widom's test particle insertion ||34|| 
proved to be sufficient. Simulation details are given in the Appendix. 

Results for vapor-liquid equilibria obtained with the new ethanol model together with ex- 
perimental data [|35l are given in Table [21 Figures [2l [3] and |4] depict the results including the 
TraPPE-UA simulation data taken from [4J . The agreement between the new model and the ex- 
perimental data is excellent, even better than that for the more complex and already very good 
TraPPE-UA ethanol model. The simulaltion results of the present ethanol model yield mean 
unsigned errors compared to experimental data [|35l in vapor pressure, saturated liquid density 
and heat of vaporization of 3.7, 0.3, and 0.9 %, respectively, in the temperature range 270 to 
490 K, which is about 55 to 95 % of the critical temperature. Following the procedure suggested 
by Lotfi et al. [f36ll . the critical temperature, density and pressure were determined. The re- 
sults compare favorably to experimental data (numbers in parenthesis): Tc =514.12 (516.25) K, 
p^. =5.94 (5.99) mol/1 and Pc =5.89 (6.38) MPa. The TraPPE-UA molecular model for ethanol 
yields mean unsigned errors in vapor pressure, saturated liquid density and heat of vaporization 
of 8.4, 1.4, and 1.6 %, respectively. Critical temperature and density are quoted with Tc =514 
K and p,. =6.1 mol/1 ffl. 
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2.3 Solutes 



Numerous Lennard- Jones based potential models for the solutes studied in the present work are 
available from the literature, e.g. for methane [I51[IIl|2ll, (3711 -(391, for nitrogen JSll-lIll, for 
oxygen (Ml l39l H IS |50l and for carbon dioxide [|30l[37l[38ll47ll.ll5ll-[|55l. 

Here, models of the symmetric two-center Lennard- Jones plus pointquadrupole (2CLJQ) 
type [|2]| were choosen, since that model class is able to describe the intermolecular interac- 
tions of all four solute molecules consistently. Furthermore, it has succesfuUy been applied to 
mixtures (SKU. The potential parameters were adjusted in previous work of our group [|2l[3]| 
exclusively to experimental pure component vapor-liquid equilibria, i.e. bubble density and va- 
por pressure. The complete parameter set of the solute models is listed in Tabled A detailed 
definition of the symmetric elongated quadrupole with momentum Q is given in |l3l. The mean 
unsigned error of the vapor pressure, saturated liquid density and heat of vaporization is smaller 
than 1, 3 and 3 %, respecitvely, for all four solute models practically over the whole vapor- liquid 
equilibrium temperature range. 

2.4 Mixtures 

In molecular simulations of binary mixtures with pairwise additive potentials, three different 
interactions are present: two between molecules of the same component which are fully defined 
by the pure substance models and the unlike interaction between molecules of different kind. 
The unlike polar interactions were determined in a straightforward manner without using bi- 
nary parameters. To define the unlike Lennard- Jones interactions between solvent and solutes 
methane, nitrogen, oxygen and carbon dioxide, the modified Lorentz-Berthelot combining rule 
was used 
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Oaa+Obb ... 
<^ab = 2 ' 

^ab = ^ ■ V^acfyb , (5) 



where ^ is the binary interaction parameter which accounts for mixture effects. 

In the "predictive mode", the Henry's law constants were computed setting ^ = 1, cf. Equa- 
tion ©. 

Simulations of many binary mixtures |I3] [H have shown that the Lorentz-Berthelot combin- 
ing rule (^=1) is too crude to predict the unlike Lennard- Jones energy parameter in an accurate 
way. Therefore, in the "adjusted mode", ^ was fitted to one experimental Henry's law constant. 
Compared to using polarizable models, the adjustment of the state independent parameter ^ is a 
simple procedure [8|. The results for the predictive and the adjusted mode calculating Henry's 
law constants are discussed in Section IH 

Other combining rules were suggested e.g. by Kohler [[591 , Hudson and McCoubrey [60 J , 
Fender and Halsey Hiza et al. [|62l-[l6l, Smith and Kong [IMlill, Waldman and Hagler 
W7\ . Al-Matar and Rockstraw [|68][ . They were mainly developed on the basis of noble gas data. 
The unlike Lennard- Jones energy parameters of binary noble gas mixtures determined from the 
above mentioned combining rules are mostly smaller than those given by the geometric mean 
of the Lorentz-Berthelot combining rule. But as shown in Section |4l the unlike solvent-solute 
interaction determined from the Lorentz-Berthelot combining rule has to be increased for three 
out of the four studied solutes > 1) to yield quantitative results. Therefore, the alternative 
combining rules will not give better predictions and were not studied further. 
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3 Henry's law constants 



Several approaches have been proposed in the literature ll56l[57ll to obtain Henry's law constants 
on the basis of molecular models. Henry's law constants are related to the residual chemical 
potential of the solute i at infinite dilution fji^ [581 

Hi = pkBTexp{ijr/ikBT)), (6) 

where is the Boltzmann constant, T the temperature, p the density of the solvent. 

In order to evaluate molecular dynamics simulations applying Widom's test particle 
method Il34l are sufficient here. This is due to the fact, that the solute molecules are all smaller 
than ethanol molecules and so acceptable statistics can be achieved. Therefore, test particles 
representing the solute i were inserted after each time step at random positions into the liquid 
solvent and the potential energy between the solute test particle and all solvent molecules \|/, 
was calculated within the cut-off radius 

= -kBT\n{Vcxp{-^\fi/{kBT)))/{V), (7) 

where V is the volume and the brackets represent the NpT ensemble average. 

The residual chemical potential at infinite dilution and hence, the Henry's law constant is 
directly attributed to the unlike solvent-solute interaction and indirectly to the solvent-solvent 
interaction which yields the configurations of the solvent molecules. In these configurations, the 
solvent test particles are inserted. The mole fraction of the solute in the solvent is exactly zero, 
as required for infinite dilution, since the test particles are instantly removed after the potential 
energy calculation. Simulations were performed at specified temperature and the according 
vapor pressure of pure ethanol. 
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4 Results 



Henry's law constants of the four solutes methan, nitrogen, oxygen and carbon dioxide in 
ethanol were determined from simulation in the predictive mode at temperatures between 273 
to 498 K with an increment of 25 K. These simulation results are given in Table IH 

For assessing the simulation results reliable experimental data are necessary. For each of 
the four studied systems, the available experimental data were grouped into two categories: 
data with low scattering that is confirmed by data sets of other authors on one hand and the rest 
of the data on the other hand. Data of the first category are called "confirmed data" here for 
brevity. Experimental data and simulation results are depicted in Figures [5] to [81 The numerical 
values for the simulation results both predictive and adjusted mode are given in Tables |4] and [51 

The examination of the experimental data at 300 K shows, that the Henry's law constants 
differ from about 150 (carbon dioxide) up to 3000 bar (nitrogen), i.e. more than one order 
of magnitude. Compared to this, the prediction of the gas solubility from pure component data 
alone are good, especially when taking into account that this property is dominated by the unlike 
interaction in mixtures. In the following, first the results of the predictive mode are discussed 
for all four solutes. After this follows the discussion for the results of the adjusted mode. 

For methane, cf. Figure [51 the predictions agree well with the experimental data. Deviations 
are about 20 %. It is remarkable that the temperature dependence of Henry's law constant of 
methane in ethanol, which shows a maximum at about 350 K, is predicted correctly. Since 
the slope of the Henry's law constant as a function of temperature H{T) is related to the heat 
of solution, this means that also that caloric property is predicted well. Also for nitrogen, the 
agreement between the predictions and the experimental data is typically within 20 %, cf. Figure 
[6l Even better results are achieved for oxygen, cf. Figure [71 where deviations are mostly within 
10 %. Only for carbon dioxide, the predicted results are worse, cf. Figure[Hl The small Henry's 
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law constants of carbon dioxide in ethanol are predicted too high by up to about 70 %. But even 
for that system, the predictions from molecular simulation lie within the range of the largely 
scattering experimental data. In all cases the predictions from the simulation for the slope of 
H{T) are excellent. Considering the simplicity and empirical nature of the Lorentz-Berthelot 
combining rule used here, this agreement of the predictions of the Henry's law constants with 
the experimental data is astonishingly good. 

As explained above, combining rules have only a weak theoretical basis and work best when 
used in a form that allows an adjustment to experimental mixture data. For that adjustment one 
single reliable mixture data point is often sufficient, even when a wide range of states of the 
mixture needs to be modelled. All this can be demonstrated very well also for the case study of 
the present work. For all four solutes, reliable experimental data on the Henry's law constants 
are available at around room temperature. Data at other temperatures are usually questionable, 
except perhaps for nitrogen. Therefore, in the present work for all four solutes, the binary 
interaction parameter ^ of the modified Lorentz-Berthelot combining rule, cf. Equation ([5]), 
was adjusted to experimental data of the Henry's law constants at 298 K (methane [[69ll , nitrogen 
fTOl , oxygen [[Tl;] and carbon dioxide (721). The adjusted binary interaction parameters are 
^ = 1.0403 (methane), 1.0449 (nitrogen), 0.9820 (oxygene) and 1.0790 (carbon dioxide) and 
were used for predictions of the Henry's law constants at other temperatures. The results given 
in Table [5] and shown in Figures [5] to [8] are what we hold to be presently the most reliable and 
available information on H{T) of the studied solutes in ethanol. The predictions are confirmed 
in the only case were the experimental data base allows such a test, i.e. for nitrogen (Figure[6l). 

Recently recommended experimental data were published on the Internet ['961. In Figures 
[5] to [HI these data were added together with their uncertainties which are typically about 5 %. 
All simulation results of the descriptive mode agree with the recommended experimental re- 
sults within the uncertainties. Hence, the estimated deviations of the Henry's law constants 
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determined in the predictive mode are validated by these recommended data. 

Additionally to the calculation of the Henry's law constants, the vapor-liquid equilibria of 
the mixture carbon dioxide and ethanol at 323 K were simulated with the Grand Equilibrium 
method ||971 using the adjusted binary interaction parameter = 1.0790). The simulation re- 
sults and the experimental data are given in Figured 

The excellent agreement of the decriptive mode results for both the Henry's law constants 
and the vapor-liquid equilibria over the whole composition range confirm that reliable mixture 
properties can be obtained over a wide range of state points with simple molecular models when 
the unlike interactions are adjusted to a small experimental data basis. 

5 Conclusion 

A new rigid molecular model for ethanol was developed modelling dispersion and repulsion 
with three Lennard- Jones sites for the methyl, methylene and hydroxyl group and modelling the 
polarity and hydrogen bonding with three point charges located at the methylene and hydroxyl 
Lennard- Jones sites, as well as at the nucleus of the hydroxyl hydrogen. The parameters were 
optimized to give accurate pure component saturated liquid density and vapor pressure. This 
potential model yields mean unsigned errors in vapor pressure, saturated liquid density and heat 
of vaporization of 3.7, 0.3 and 0.9 %, respectively, in the temperature range from 270 to 490 
K. This new molecular model for ethanol was applied together with two-center Lennard- Jones 
plus pointquadrupole (2CLJQ) models [|2| for methane, nitrogen, oxygen and carbon dioxide to 
predict Henry's law constants. 

Using the Lorentz-Berthelot comining rule, simulations were performed to predict Henry's 
law constants of the four solutes in ethanol from 273 to 498 K with an increment of 25 K. For 
methane, nitrogen and oxygen, favorable predictions are achieved with deviations below 20 %. 



14 



For carbon dioxide, for which the Henry's law constants are small, deviations of up to 70 % are 
observed. 

As an alternative to predict from pure component data alone, the modified Lorentz-Berthelot 
combining rule was studied. The binary interaction parameter ^ was adjusted to one experi- 
mental Henry's law constant at 298 K, a temperature for which reliable data on the Henry's 
law constants in ethanol are available for all studied solutes. The models with the adjusted bi- 
nary interaction parameter give what we believe is the presently most reliable information on 
H{T) in the studied systems. The adequacy of the mixture models are supported by the vapor- 
liquid equilibrium simulation results for the mixture carbon dioxide and ethanol over the whole 
composition range. 

The methods applied in this work to calculate Henry's law constants are limited to solute 
molecules of smaller size than the solvent molecules due to Widom's test particle method. In 
future work, the feasibility with strongly dipolar solute molecules will be examined. 
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7 Appendix, Simulation Details 

In all simulations A^^ = 864 molecules were used. The calculation of the vapor-liquid equilibra of 
pure ethanol was done with the NpT+test particle method [|28l - [|30l . In order to determine the 
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chemical potential and pure ethanol properties in the liquid phase, the gradual insertion method 
||3T| - [[33l was applied. In such a simulation run, the liquid is equilibrated in the NpT ensem- 
ble sufficiently over 10 000 cycles without fluctuating molecules and adjusting the maximum 
displacement of translation, rotation and volume to yield acceptance rates of 50 %. After that, 
5 000 Monte Carlo cycles with fluctuating molecules are performed to equilibrate the weights of 
the transition states. The production phase was performed over 50 000 Monte Carlo cycles. The 
nonentropic properties and the chemical potential in the vapor phase were determined applying 
standard Monte Carlo simulation and Widom's test particle insertion ll34l . A vapor simulation 
consists of a sufficiently long equilibration phase and 200 000 Monte Carlo cycles where test 
particles are inserted after each cycle. 

In order to evaluate the Henry's law constants, molecular dynamics simulations have been 
performed. After a long equilibration of 200 000 time steps, 250 000 production time steps 
were carried out inserting AN test particles after each step. 

The Lennard- Jones long range interactions beyond the cut-off radius were corrected em- 
ploying angle averaging proposed by Lustig [17311 . The coulombic interactions were corrected 
using the reaction field method JTH. The cutoff radius was at least 17.5 A. The quadrupolar 
interaction needs no long range correction, as its contribution disappears by orientational aver- 
aging. For the simulation of pure ethanol vapor, the diameter of the cutoff sphere was chosen 
close to half simulation box length due to faster equilibration and better sampling. 
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Table 1: Lennard- Jones, point charge and geometry parameters of the present ethanol model, 
cf. Equation ([T]) and Figure[Il electronic charge e = 1.602177 ■ 10^^^ C. 



Site 


<3aa 

A 


K 


qia 

e 


ScH3 


3.6072 


120.15 





ScH2 


3.4612 


86.291 


0.25560 


SoH 


3.1496 


85.053 


-0.69711 


Sh 








0.44151 



hi 


hi 




Yi 


Y2 


A 


A 


A 


deg 


deg 


1.98420 


1.71581 


0.95053 


90.950 


106.368 



Table 2: Vapor-liquid equilibrium of pure ethanol: simulation results compared to experimental 
data ll35l for vapor pressure, saturated densities and heat of vaporization. The numbers in 
parenthesis indicate the statistical uncertainty in last digit. 



T 
K 


Psim 

MPa 


Pexp 

MPa 


Psim 

mol/1 


Pexp 

mol/1 


p" 

r Sim 

mol/1 


o" 

r exp 

mol/1 


Ah^- 

Sim 

kJ/mol 


Ah^ 

'-"'exp 

kJ/mol 


270 


0.00125(4) 


0.00129 


17.68(1) 


17.64 


0.0008(4) 


0.00056 


45.03(4) 


44.78 


314 


0.01951(7) 


0.01868 


16.83(1) 


16.76 


0.0077(4) 


0.00710 


42.03(4) 


42.33 


358 


0.1239 (5) 


0.12977 


15.84(1) 


15.78 


0.044 (1) 


0.04514 


38.32(4) 


38.80 


402 


0.526 (3) 


0.55451 


14.66(1) 


14.65 


0.180 (1) 


0.18623 


33.40(6) 


33.99 


446 


1.60 (3) 


1.69564 


13.16(2) 


13.24 


0.56 (1) 


0.58896 


27.1 (1) 


27.48 


490 


3.88 (8) 


4.08339 


11.10(4) 


11.14 


1.5 (2) 


1.68190 


19. (1) 


17.96 
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Table 3: Parameter set of the 2CLJQ potential for the solutes, taken from 



Solute 


a 




L 


Q 




A 


K 


A 


DA 


CH4 


3.7281 


148.55 







N2 


3.3211 


34.897 


1.0464 


-1.4397 


O2 


3.1062 


43.183 


0.9699 


-0.8081 


CO2 


2.9847 


133.22 


2.4176 


-3.7938 



Table 4: Henry's law constants of four solutes in ethanol for the predictive mode = 1). The 
numbers in parenthesis indicate the statistical uncertainty in the last digits. 



T 


HcH4 


Hm 


H02 


Hco2 


K 


bar 


bar 


bar 


bar 


273 


923 (16) 


3647 (56) 


1513 (19) 


221.1 (38) 


298 


997 (14) 


3381 (43) 


1542 (15) 


302.4 (39) 


323 


1076 (12) 


3177 (32) 


1593 (12) 


390.2 (39) 


348 


1069.3 (79) 


2810 (20) 


1540.0 (84) 


456.5 (30) 


373 


1070.9 (67) 


2536 (16) 


1491.8 (73) 


515.3 (29) 


398 


1007.8 (57) 


2171 (12) 


1374.7 (64) 


540.7 (28) 


423 


929.9 (58) 


1835 (12) 


1236.9 (73) 


547.2 (31) 


448 


817.8(39) 


1483.0 (78) 


1068.1 (46) 


524.1 (23) 


473 


703.2 (44) 


1172.4 (84) 


896.3 (54) 


485.5 (27) 


498 


552.1 (43) 


837.4 (77) 


688.4 (55) 


411.9 (28) 
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Table 5: Henry's law constants of four solutes in ethanol for the adjusted mode: methane 
^=1.0403, nitrogen ^=1.0449, oxygen ^=0.9802 and carbon dioxide ^=1.0790. The numbers 
in parenthesis indicate the statistical uncertainty in the last digits. 



T 


Hem 


^N2 


Ho2 




K 


bar 


bar 


bar 


bar 


273 


686 (12) 


2867 (43) 


1691 (17) 


101.7 (17) 


298 


797 (11) 


2827 (34) 


1765 (16) 


158.1 (20) 


323 


848.4 (89) 


2614 (26) 


1763 (14) 


213.8 (20) 


348 


877.1 (82) 


2390 (21) 


1691 (11) 


269.1 (22) 


373 


878.0 (59) 


2148 (14) 


1601.0(99) 


319.5(19) 


398 


865.0 (51) 


1942 (11) 


1483.4(69) 


360.3 (19) 


423 


809.4 (43) 


1636 (11) 


1306.1 (70) 


381.2(20) 


448 


733.7 (38) 


1343.2 (69) 


1117.9(49) 


383.6(15) 


473 


634.0 (38) 


1078.2 (75) 


931.6 (57) 


372.7 (20) 


498 


507.9 (38) 


783.3 (70) 


705.0 (57) 


332.8 (20) 
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